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We  have  derived  an  analytic  model  describing  the  interior  temperature  difference  as  a  function  of  the 
load  current  of  a  thermoelectric  generator  (TEG);  we  have  also  proposed  a  method  to  extract  the  intrinsic 
and  extrinsic  Seebeck  coefficients  and  resistances  of  TEG  using  experimental  current-voltage  curves.  The 
decrement  of  internal  temperature  difference  is  almost  linearly  depending  on  load  current  of  the  TEG. 
From  the  experimental  results,  using  a  thermoelectric  (TE)  module  with  a  thermal  conductance  of 
3.52  W/K  and  a  parasitic  thermal  conductance  of  50  W/K,  the  effective  internal  electrical  resistance 
was  increased  by  approximately  5%,  but  the  effective  Seebeck  coefficient  was  attenuated  by  approxi¬ 
mately  13%,  as  compared  to  the  intrinsic  parameters.  The  relationship  between  the  output  power  reduc¬ 
tion  factor  and  limited  thermal  conductance  of  TEG  packaging  was  also  derived.  Approximately  25%  of 
the  maximum  output  power  is  lost  because  of  the  parasitic  thermal  resistance  of  the  TE  module  used 
in  the  experiment. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  recycling  of  waste  heat  to  electrical  energy  using  a  thermo¬ 
electric  generator  (TEG)  has  recently  attracted  attention  in  areas 
where  a  considerable  amount  of  energy  is  wasted,  such  as  vehicu¬ 
lar  applications  [1-8].  The  TEG  is  a  solid-state  energy  converter 
that  generates  electric  power  from  the  temperature  difference 
between  its  two  sides.  The  internal  temperature  difference  and 
its  dependence  on  the  load  current  have  not  been  characterized 
in  analytic  models  [7,8].  The  analysis  and  modeling  of  the  effective 
temperature  difference  and  electrical  properties  of  the  TEG  are 
very  important  issues  in  the  design  of  an  optimum  waste  heat 
recovery  system  [9].  The  prior  arts  of  one-dimensional  models 
can  be  applied  to  the  study  of  internal  temperature  of  TEG 
[10-13].  However,  as  long  as  the  author  knows,  an  extraction 
method  of  model  parameters  of  thermoelectric  generator  while 
load  current  is  flowing  has  not  been  reported. 

In  a  thermoelectric  (TE)  module,  p-type  and  n-type  thermoelec¬ 
tric  pellets  sandwiched  between  ceramic  plates  are  electrically 
connected  in  series  with  a  Cu  film,  as  shown  in  Fig.  1  [10-16]. 
The  effective  thermal  resistance  of  the  ceramic  plate  and  Cu  inter¬ 
connection  layer  on  both  sides  of  the  TE  module  (0H  and  6C)  and  the 
thermal  resistance  of  the  pellets  (0m)  are  connected  in  series  and 
form  a  thermal  resistance  network  [3,17,18].  Generally,  a  TE  mod¬ 
ule  is  a  symmetric  structure,  designed  to  have  an  identical  thermal 
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resistance  between  both  hot  or  cold  plates  and  pellets.  The  thermal 
and  electrical  characteristics  of  the  TEG  are  represented  by  the 
SPICE  equivalent  model  shown  in  Fig.  2  [12,14-16].  Thermal  resis¬ 
tances  9h  or  0C  include  parasitic  thermal  resistance  caused  by  the 
ceramic  plate,  Cu  layer,  and  thermal  interface  material  (TIM),  and 
0m  is  the  effective  thermal  resistance  formed  by  a  parallel  combina¬ 
tion  of  all  of  the  pellets. 

TEG  can  be  modeled  as  a  Thevenin  voltage  source  proportional  to 
the  Seebeck  coefficient  with  an  intrinsic  temperature  difference 
between  its  hot  and  cold  sides,  am(Te  -  Ta),  and  an  intrinsic  resis¬ 
tance  of  Rmi  as  shown  in  Fig.  2b  [  1 4-1 6].  The  external  generated  volt¬ 
age-current  relationship  is  modeled  as  a  Thevenin  voltage  source 
proportional  to  the  effective  Seebeck  coefficient,  with  a  temperature 
difference  between  the  hot  and  cold  plates,  ocTeg(Th  -  Tc),  and  an 
effective  resistance  of  RTEG ,  as  shown  in  Fig.  2c.  The  generated  output 
voltage-current  relationship  is  represented  by  [12,14-16], 

V  —  0Cteg(Th  —  Tc)  —  RtecI  =  MtegAT  —  RtegT  (1) 

In  conventional  analysis,  the  effective  Seebeck  coefficient  (ocTEg) 
is  extracted  from  the  measured  open  circuit  voltage  divided  by  the 
external  temperature  difference  (TH  -  Tc),  and  RTEG  is  extracted 
from  the  slope  of  the  I-V  curve.  If  the  thermal  resistance  of  the 
ceramic  plates  and  Cu  interconnect  is  negligible,  the  temperature 
of  the  pellets  is  identical  to  the  external  temperature,  but  in  the 
TE  module,  the  internal  temperature  difference  of  the  pellets 
(Te  -  Ta)  is  lower  than  the  external  temperature  difference 
(Th  -  Tc).  Because  of  the  parasitic  thermal  resistance,  the  tempera¬ 
ture  difference  diminished  even  under  the  open  circuit  condition  at 
the  output  node  of  the  TEG  [17,18]. 


0306-2619/$  -  see  front  matter  ©  2012  Elsevier  Ltd.  All  rights  reserved. 
http://dx.doi.org/!  0.1 01 6/j.apenergy.201 2.09.006 


S.  Kim/ Applied  Energy  102  (2013)  1458-1463 


1459 


Hot  side 


TE  Pellets 


Cold  side 

Fig.  1.  Structure  and  thermal  resistor  network  of  thermoelectric  module. 
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2.  Modeling  a  thermoelectric  generator 

The  rate  of  supply  heat  Qh  and  removal  heat  Qc  can  be  esti¬ 
mated  at  the  hot  and  cold  junction  respectively.  The  heat  flow 
equations  of  each  side  of  the  hot  and  cold  plates  are  [11-16], 

t2r 

Qh  =  km(Te  -  To)  +  XmlTe - ^  =  kc(TH  -  Te),  (2) 

Qc  =  l<m(Te  -  Ta)  +  amITa  +  ^  =  kc(Ta  -  Tc).  (3) 

Here,  (Te  -  Ta)  is  the  temperature  difference  of  the  internal  plate 
of  the  p-  and  n-type  pellets,  and  (TH  -  Tc )  is  the  temperature  differ¬ 
ence  of  the  external  plate  of  the  TEG.  kc  and  km  are  the  effective 
thermal  conductances  of  the  parasitic  components  and  pellets, 
respectively,  and  are  the  inverse  values  of  thermal  resistances  9C 
and  0m,  respectively. 

The  device  parameters  of  TEG  linearly  depend  on  area  (A), 
length  (L)  and  the  number  of  pellets  (N)  [13].  These  parameters 
can  be  calculated  as 

.  2AM  .  2 AM  _  2 NL  , 

km  —  K  ^  5  kc  —  K(Cu+Ceramic)  ^  5  —  P  ^  1  and  0Lm  —  (2 1V)S,  (4) 

where  s,  k  and  p  are  the  Seebeck  coefficient,  thermal  conductivity 
and  electrical  resistivity  of  the  pellets,  respectively.  K(Cu+ceramic )  is 
effective  thermal  conductivity  of  parasitic  components  composed 
of  the  ceramic  plate  and  Cu  interconnects. 

We  assume  a  symmetrical  structure  for  the  hot  and  cold  plates, 
so  the  parasitic  thermal  conductance  of  the  ceramics  and  Cu  plates 
are  equivalent  for  the  hot  and  cold  plates  of  the  TEG. 

By  subtracting  Eq.  (3)  from  Eq.  (2), 
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Fig.  2.  Electro-thermo  equivalent  circuit  model  of  TEG.  (a)  Electric  circuit  model  of 
thermal  behavior,  (b)  Equivalent  source  model  of  TEG  with  internal  temperature 
difference  and  (c)  Thevenin  equivalent  source  model  of  TEG  with  external 
temperature  difference. 


a J(Te  -  Ta)  -  I2Rm  =  kc(TH  +  Tc)  -  kc(Te  +  Ta).  (5) 

By  adding  Eq.  (2)  and  (3), 

2 km(Te  -  Ta)  +  2ocmI(Te  +  Ta)  =  kc(TH  -  Tc)  -  kc(Je  -  Tfl),  (6) 

am  and  Rm  are  the  intrinsic  Seebeck  coefficient  and  electrical  resis¬ 
tance,  respectively,  of  the  TEG  used  to  model  the  electrical  charac¬ 
teristics  of  the  TEG. 

From  Eqs.  (5)  and  (6),  we  obtain  an  expression  for  the  temper¬ 
ature  difference  between  the  internal  hot  and  cold  plates  of  the 
TEG  (AT/nt), 


It  is  quite  difficult  to  measure  the  internal  temperature  difference 
between  the  ends  of  the  pellets  while  the  load  current  is  flowing.  Un¬ 
der  the  condition  of  a  nonzero  load  current,  the  internal  temperature 
difference  drops,  even  though  the  external  temperature  remains 
constant.  In  prior  works  [17,18],  the  actual  internal  temperature  dif¬ 
ference  is  calculated  by  including  the  parasitic  thermal  resistance  of 
the  TE  module  under  the  open  circuit  condition,  and  the  effective 
Seebeck  coefficient  was  obtained  from  the  experimental  results. 
However,  while  the  load  current  is  flowing,  the  internal  temperature 
difference  of  the  TEG  drops,  which  diminishes  its  power  generation 
efficiency  [9].  There  have  been  no  prior  reports  regarding  the  ana¬ 
lytic  modeling  of  the  internal  temperature  difference  as  a  function 
of  the  load  current  or  on  the  influence  of  the  accuracy  of  the  ex¬ 
tracted  electrical  parameters  and  power  generation  efficiency. 

In  this  paper,  we  derive  an  equation  describing  the  interior  tem¬ 
perature  difference  as  a  function  of  the  output  load  current  of  the 
TEG  and  propose  a  method  to  extract  the  effective  Seebeck  coeffi¬ 
cient  and  electrical  resistance  of  the  TEG.  Theoretical  analysis  and 
parameter  extraction  methods  from  the  measured  data  and  valida¬ 
tion  of  the  approach  using  SPICE  circuit  simulation  results  will  be 
explained  in  the  following  sections. 


C Te  ~  Ta)  = 


kc 


2 kmkc  +  k2c  -  oc^I2 
=  AT/nt . 


kc(TH  -  Tc)  - 


~~kT 


-  ocmI{TH  +  Tc) 


(7) 


The  Eq.  (7)  is  obtained  from  the  heat  flow  equations,  Eqs.  (2)  and 
(3)  at  hot  and  cold  plate,  respectively.  Whereas,  we  assume  that  the 
external  temperature  difference  between  hot  and  cold  plate, 
(Th  -  Tc),  is  kept  constant  during  operation.  The  Eq.  (7)  indicates 
that  there  are  three  terms  influencing  the  internal  temperature  dif¬ 
ference.  The  first  factor  is  the  voltage  divide  ratio  by  thermal  resis¬ 
tive  network.  The  second  and  third  terms  are  representing  wattage 
caused  by  Joule  heating,  and  wattage  due  to  Peltier  effect  respec¬ 
tively,  while  load  current  is  flowing. 

When  the  load  current  is  zero,  we  obtain  A TInt  under  the  open 
circuit  condition. 


AT/nt  |/^o 


=  ATj0  = 


kc(TH  -  Tc) 

2km  +  kc 


(8) 


This  expression  is  the  same  as  the  internal  temperature  differ¬ 
ence  calculated  using  the  thermal  resistor  network  of  a  single  TE 
couple  [17,18]. 
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By  differentiating  A TInt  with  respect  to  load  current  /,  we  obtain 
the  slope  of  the  change  of  the  temperature  difference,  0(  A T/nt)/0J,  as 
follows: 

<9(AT/nt)  ^  Oim(jH  +  Tc) 

~d^=  2km  +  kc  •  [  ) 

The  derivation  of  Eq.  (9)  is  explained  in  Appendix  A.  The  phys¬ 
ical  meaning  of  Eq.  (9)  is  that  internal  temperature  difference  of 
the  TEG  can  be  approximated  by  a  linear  function  of  load  current, 
and  Peltier  effect  is  the  major  factor  of  internal  temperature  drop 
caused  by  the  current  flowing  through  a  TEG.  The  slope  of  change, 
0(ATjnt)/0J,  is  dependent  with  size  of  TE  module,  because  the  value 
of  TE  module  parameters,  km,  kc  and  am,  are  linear  function  of  area, 
length  and  the  number  of  pellets  [13]. 

Because  the  slope  of  the  change  of  the  temperature  difference  is 
almost  independent  of  the  load  current,  as  shown  in  Eq.  (9),  we  can 
model  the  linear  dependence  of  the  interior  temperature  difference 
on  the  load  current. 


J<c(Th-Tc)  qm(TH  +  Tc)/ 
2km  +  kc  2km  +  kc 


(10) 


The  first  term  in  Eq.  (9)  is  the  temperature  attenuation  due  to 
the  thermal  resistor  network  comprising  kc  and  kmi  and  the  second 
term  is  the  temperature  drop  caused  by  the  Peltier  effect.  From  the 
electrical  model  of  the  TEG,  we  obtain  the  relationship  between  the 
terminal  voltage  and  load  current,  as  follows: 

V  =  ocmATlnt  -  RJ  -  <xmATi0  -  ( Rm  +  ^^±^1 

=  2^TFc(Th~Tc)~RtecI-  (11) 


where  Rm  is  the  internal  resistance  and  Rteg  is  the  effective  terminal 
electrical  resistance  of  the  thermoelectric  generator. 

From  Eqs.  (1)  and  (11),  we  can  determine  the  relationship  be¬ 
tween  the  intrinsic  and  extrinsic  parameters  of  the  TEG, 


Mteg  = 


kc 


2km  +  kc 


(12) 


and  Rteg 


<(Th  +  Tc)Y 

2  km  +  kc  J 


(13) 


3.  Extraction  of  TEG  parameters 

A  TEG  module  (TMH400302055,  manufactured  by  Wise  Life 
Technologies)  with  199  TE  couples  was  used  in  the  experiment. 
The  geometric  features  and  thermal  and  electrical  parameters  of 
the  TE  module  used  in  the  experiments  are  listed  in  the  Refs. 
[3,18].  The  calculated  thermal  conductance  (km)  of  all  398  pellets 
of  TMH400302055,  including  the  TIM  between  the  Cu  interconnect¬ 
ing  layer,  is  approximately  3.52  W/K,  and  the  thermal  conductance 
of  the  ceramic  plate  and  Cu  interconnects  (/<c)  is  approximately 
50.0  W/K. 

Fig.  3  shows  the  experimental  curves  of  the  output  voltage 
versus  the  current,  which  are  fit  well  to  straight  lines  for  several 
measurement  condition  of  ( TH  -  Tc),  as  reported  in  prior  works 
[1,3,10-17]. 

We  have  no  direct  approach  for  measuring  the  interior  temper¬ 
ature  difference  A TInt  of  the  TEG  during  power  generation.  We 
therefore  propose  a  procedure  for  the  extraction  of  the  interior 
temperature  difference  as  well  as  the  electrical  parameters  using 
the  measured  electrical  data,  as  follows: 

(1)  Extraction  of  thermal  conductance,  km  and  km. 

(2)  Calculation  of  A Ti0  using  Eq.  (8). 


Current  (A) 

(a) 


Current  (A) 


(b) 


Current  (A) 

(C) 


Fig.  3.  Measured  voltage-versus-current  characteristics  of  a  Bi2Te3  thermoelectric 
module  for  various  hot  side  temperatures  ( TH )  and  temperature  differences 
(AT  =Th-  Tc)  between  hot  and  cold  plates,  (a)  For  TH  =  320  K,  (b)  for  TH  =  350  K, 
and  (c)  for  T„  =  370  K. 


(3)  Extraction  of  ocTeg  and  RTEG  from  the  V-axis  intersection  (Voc) 
and  slope  of  the  curve,  respectively.  From  the  measured  I-V 
curve,  ocTEG  =  V0J(TH  -  Tc)  and  RTEG  is  the  slope  of  the  line  of 
best  fit  of  the  measured  data. 

(4)  Calculation  of  am,  using  Eq.  (12) 

(5)  Calculation  of  Rm ,  using  Eq.  (13) 

(6)  Calculation  of  A TInt,  using  Eq.  (7)  or  Eq.  (10) 
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We  only  need  two  thermal  conductance  parameters,  kc  and  km, 
for  the  analysis  and  extraction  of  the  electrical  parameters  from 
the  measured  data. 


4.  Experimental  results  and  discussions 


We  calculated  the  internal  temperature  difference,  Seebeck 
coefficient,  and  internal  electrical  resistance  according  to  the  pro¬ 
posed  procedure,  using  the  experimental  I-V  curves.  The  calculated 
internal  temperature  difference,  as  a  function  of  load  current  using 
Eqs.  (7)  and  (10),  shows  a  good  linear  fit.  We  found  the  slope  of  the 
change  of  temperature  difference  A TInt  with  respect  to  the  load 
current  to  be  approximately  -0.88  °C/A,  which  coincided  with 
the  value  calculated  using  Eq.  (10),  as  shown  in  Fig.  4.  The  experi¬ 
mental  I-V  curves  for  hot  side  temperatures  of  330  K  and  340  K  are 
not  shown  in  this  paper  for  simplicity,  because  the  trend  of  the 
experimental  curve  was  similar.  The  thermal  resistance  of  the  con¬ 
tact  material  comprising  the  ceramic  plate,  Cu  interconnect,  and 
solder  between  the  TE  pellets  and  Cu,  degrade  the  performance 
of  the  TEG  [17,18].  Fig.  5  shows  the  extracted  values  of  the  intrinsic 
and  external  Seebeck  coefficients  as  functions  of  temperature.  The 
external  Seebeck  coefficient  is  approximately  87%  of  the  intrinsic 
one,  a  degradation  ratio  that  is  consistent  with  the  results  of  prior 
work  using  the  same  type  of  TEG  [17,18]. 

Fig.  6  shows  the  extracted  electrical  resistance  obtained  from 
the  slopes  of  the  measured  I-V  curves,  RTEG,  and  the  calculated 
internal  resistance  value,  Rm,  using  Eq.  (13).  The  effective  external 
resistance  value  is  increased  by  approximately  5%  because  of  a 
drop  in  the  internal  temperature  difference  in  the  temperature 
range  320-370  K  in  the  experiment  when  kc  is  approximately 
50W/K  and  km  is  approximately  3.52  W/K.  Experimental  results 
from  Refs.  [14-19]  shows  that  temperature  dependence  of  Seebeck 
coefficient  and  electrical  resistivity  of  bulk  thermoelectric  material 
is  nonlinear.  The  electrical  conductivity  and  Seebeck  coefficient  are 
tightly  related  to  the  carrier  concentration,  and  carrier  mobility 
caused  by  various  scattering  effects  which  strongly  depends  on 
temperature  [10,19]. 

Fig.  7  shows  the  measured  data,  the  SPICE-simulated  data  using 
the  equivalent  circuit  shown  in  Fig.  2b,  and  the  calculated  data 
from  Eq.  (10).  The  calculated  and  simulated  data  using  the  pro¬ 
posed  model  and  extraction  procedure  coincide  quite  well  with 
the  measured  data. 

Because  the  internal  electrical  resistance  of  the  TEG  in  the  SPICE 
model  includes  thermal  behavior,  it  is  necessary  to  use  Rm  instead  of 
Rteg •  Rteg  implicitly  includes  the  load  current  effect  on  the  drop  in 
internal  temperature  difference.  For  this  reason,  in  the  discrete  mod¬ 
el  that  does  not  consider  the  cross-coupled  action  of  the  electro¬ 
thermal  behavior,  it  is  necessary  to  use  RTEC  as  an  internal  source 
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Fig.  4.  Calculated  internal  temperature  difference  using  Eq.  (6)  (markers)  and  Eq. 
(9)  (lines). 
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Fig.  5.  Extracted  Seebeck  coefficients  from  the  measured  I-V  curves. 
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Fig.  6.  Extracted  electrical  resistances  from  the  measured  I-V  curves. 


resistance  value.  It  is  very  hard  to  specify  the  limitation  of  tempera¬ 
ture  difference  that  the  SPICE  model  could  be  applied  on  TEG  mod¬ 
ules.  For  larger  temperature  difference,  we  need  to  extract  SPICE 
parameters  for  an  average  temperature  of  the  TE  pellet  [12-14,16]. 

The  finite  value  of  the  thermal  conductance  between  the  TE  pel¬ 
lets  and  the  contact  plate  of  the  heat  source  causes  significant  deg¬ 
radation  in  the  output  performance  of  the  TEG  because  of  an 
increase  in  the  internal  resistance  and  attenuation  of  the  effective 
Seebeck  coefficient. 

There  may  be  nonlinear  effects,  in  a  wide  range  of  temperature 
difference  between  hot  and  cold  side,  or  at  high  hot  temperature 
operation,  due  to  the  temperature  dependence  of  material  param¬ 
eters  [11-16,19].  For  example,  the  hot  temperature  can  go  from 
373  K  to  623  K  in  vehicular  applications  of  TEG  [1-8,19,20]. 

The  basic  assumption  of  the  approach  is  that  heat  flow  equa¬ 
tions,  Eqs.  (2)  and  (3),  are  valid,  and  thermal  conductance  is  not 
changed  by  load  current.  In  order  to  apply  the  proposed  procedure 
in  a  wide  range  of  temperature  difference,  we  need  to  re-calculate 
thermal  conductance  value  of  the  operating  temperature  at  the  ini¬ 
tial  step  of  the  extraction  procedure  using  datasheets  [21].  The 
main  goal  of  this  paper  is  to  derive  an  analytic  model  describing 
the  interior  temperature  difference  as  a  function  of  the  load  cur¬ 
rent.  Providing  a  simple  method  to  extract  the  intrinsic  and  extrin¬ 
sic  Seebeck  coefficients  and  resistances  of  TEG  using  experimental 
current-voltage  curves  is  another  purpose  of  the  work.  The  inves¬ 
tigation  of  physical  reason  of  temperature  dependence  of  thermal 
conductance  at  high  temperature  is  beyond  the  scope  of  this  work. 
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Fig.  7.  Voltage-current  curves  of  measured  data  (markers),  SPICE-simulated 
(dotted  lines)  data,  and  calculated  data  using  Eq.  (10)  (solid  lines),  (a)  For 
Th  =  320  K,  (b)  for  TH  =  350  K,  and  (c)  for  TH  =  370  K. 


ZT  is  the  figure  of  merit  of  the  TEG.  The  derivation  procedure  for 
Eq.  (14)  is  described  in  Appendix  B. 

In  Eq.  (14),  the  theoretical  maximum  power  generated  by  the 
TEG  is 


«m(TH  -  Tcf 

4Rm 


(16) 


Hence,  the  performance  derating  factor  due  to  the  limited  ther¬ 
mal  conductance  of  the  plates  is 

yi  —  _ (fcc/fcm) _ 

ld  (2  +  kc/km){(2  +  kc/km)  +Zr/4} '  1  J 

For  the  given  values  of  kc  =  50  W/K  and  km  =  3.52  W/K,  the 
output  power  derating  factor,  rjd,  is  approximately  0.75,  implying 
that  the  output  power  is  reduced  by  25%  owing  to  the  limited  ther¬ 
mal  conductance  of  the  packaging  and  the  interconnecting  mate¬ 
rial  comprising  the  ceramic  plate  and  the  Cu  interconnect  layer 
of  the  TE  module.  The  derating  factor  is  a  function  of  geometry  of 
a  TEG,  because  the  thermal  conductances  of  TEG  linearly  depend 
on  the  area,  length,  and  the  number  of  TE  module. 


5.  Conclusions 

We  have  derived  an  equation  describing  the  interior  tempera¬ 
ture  difference  as  a  function  of  the  load  current  of  the  TEG  and 
have  proposed  a  method  to  extract  the  effective  Seebeck  coefficient 
and  electrical  resistance  of  the  TEG  using  experimental  I-V  curves. 
The  internal  temperature  difference  decreases  almost  linearly  with 
the  load  current  of  the  TEG  because  the  wattage  due  to  Peltier 
effect  creates  temperature  drop  across  interface  materials.  From 
the  experimental  results  using  the  TE  module  with  a  thermal  con¬ 
ductance  of  3.52  W/K  and  a  parasitic  thermal  conductance  of 
50  W/K,  the  effective  internal  electrical  resistance  increased  by 
approximately  5%,  but  the  effective  Seebeck  coefficient  was  atten¬ 
uated  by  approximately  13%,  as  compared  to  the  intrinsic  parame¬ 
ters.  The  theoretical  relationship  of  the  output  power  reduction 
factor  caused  by  internal  temperature  drops  was  also  derived. 
Approximately  25%  of  the  maximum  output  power  was  lost  be¬ 
cause  of  the  parasitic  thermal  resistance  of  the  TE  module  used 
in  the  experiment.  Parameter  extraction  and  analysis  using  the 
measured  data  and  SPICE-simulated  results  show  that  the  pro¬ 
posed  approach  is  valid  for  the  modeling  and  estimation  of  the  per¬ 
formance  of  a  TEG. 
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Maximum  output  power  is  achieved  under  the  load  matching 
condition,  when  the  load  resistance  is  equal  to  the  source  resis¬ 
tance  of  the  TEG.  Since  the  I-V  curve  is  linear,  the  maximum  output 
power  is  [1,22-25] 


P MAX  — 


4  Rti 


(14) 


By  substituting  the  relationship  of  Voc  and  RTEG  into  Eq.  (13),  we 
obtain 


Pmax 


«2m(rH-rc)2 

4  Rm 


(kc/km)2 

(2  +  kc/km){(2  +  kc/km)  +Zj/4} ' 


(15) 


Appendix  A 


Rewriting  Eq.  (7),  the  equation  of  the  internal  temperature 
difference  of  the  TEG  is 


A  TInt  =  - 


2  kmkc  +  kc  -  a  If 


c-dl3  -  el 
a 2  -  b2I2 


(kc(TH  -  Tc)  - 


MmRmI 

~~kT 


-  0CmI(TH  +  Tc) 


Variables  a,  b,  c,  d ,  and  e  in  Eq.  (Al)  are  as  follows: 


(Al) 


x 
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Table  A1 

Typical  values  of  the  TE  module  parameters  and  coefficients  used  in  the  experiment 
described  in  this  paper. 


TEG  parameters 

Variables  in  Eqs.  (A1)-(A7) 

Parameters 

Values 

Coefficients  Values 

Coefficients  Values 

kc 

50 

a 

2852 

bd 

4.1  x  10“4 

km 

3.52 

b 

4.9  x  10“3 

be 

12 

OCm 

0.07 

c 

50000 

be 

245 

Rm 

1.2 

d 

8.5  x  10“2 

ae 

7.2  x  106 

Th 

370 

e 

2520 

ab 

14 

Tc 

350 

ad 

240 

a2 

8.2  x  106 

a  —  2  kmkc  +  kc , 

(A2) 

b  = 

(A3) 

c  =  k2(TH-Tc), 

(A4) 

d  —  ccmRm, 

(A5) 

e  =  ocmkc(TH  +  Tc). 

(A6) 

Differentiating  the  temperature  difference  with  respect  to  load 
current  /  gives 


d(A T,nt)  _  bdl4  -  (3 ad  +  be)/2  +  2 bcl  -  ae 
di  b/4  -  2 ab/2  +  a2 


(A7) 


Because  the  Seebeck  coefficient  of  TEG,  am,  is  relatively  very 
small,  the  effective  value  of  parameters  b ,  be,  bd,  be,  and  d  is  negli¬ 
gible  compared  to  the  values  of  ab,  ad,  or  ae. 

As  an  example,  the  typical  values  of  the  TE  module  parameters 
and  coefficients  used  in  the  experiment  described  in  this  paper  are 
listed  in  Table  Al. 

We  can  simplify  Eq.  (A7)  as  follows: 


d(ATInt)  ^  3 adl2-hae  _  3 d/2  +  e  ^  e 

dl  =~  -2 abl2  -ha2~~  -2 bl2  -ha  =  ~a 
_  _  MmiTn  +  Tc) 

2km  +  kc 


(A8) 


Appendix  B 


By  substituting  the  relationship  of  the  open  circuit  voltage,  Voc 
and  R-teg  into  Eq.  (14),  we  obtain 


Pm  ax  = 


Ga^H-rc))2 

(ocm(TH  -  Tc))2 


(kc/2kw 


4  Rm 

,<(Th-Tc)2 


\+al(TH  +  Tc)/4Rm(2km+kc 


4R m  (2 km  +  kc)  ((2k*  +  kc)  +  ' 

Since  figure  of  merit  ZT  is  given  by  [10-16] 


ZT  = 


+  TC) 


Rmkm 

We  can  rewrite  Eq.  (Bl)  as  follows: 


a  l(JH-Tc)\, 
Pmax= - ARm - X 


(kc/km 


(2  +  kc/km){(2  +  kjkm)  +Zt/4} 


(Bl) 


(B2) 


(B3) 


In  the  case  of  a  small  ZTf  Zr/4  is  negligible  compared  to  2  +  kc/km, 
and  hence,  PMax  is  simplified  to 


P MAX 


TH^Tcf  (  kc/K  \2 

4Rm  \2  +  kc/kJ  ■ 


(B4) 
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